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Abstract. A discrete analog of the Tzitzeica equation is found in the form 
of quad-equation. Its continuous symmetry is an inhomogeneous Narita- 
Bogoyavlensky type lattice equation which defines a discretization of the 
Sawada-Kotera equation. The integrability of these discretizations is proven 
by construction of the Lax representations. 
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1 Introduction 

The celebrated Tzitzeica equation [1, 2, 3, 4] 

Hxy = e H_ e ~2H (1) 

is probably the most exotic integrable equation among those which can be 
written in such a compact form. It looks deceptively similar to the Liouville 
and the sine-Gordon equations, but the underlying spectral problem is of 
third order and this makes all properties more complicated. In particular, 
the problem of integrable discretization is rather difficult. The most natural 
pattern for a discrete version is given obviously by quad-equations, that is, 
difference equations on the square grid 

Q(h,hi,h 2 ,h 12 ) = 

where h = h(ni,ri2) and the subscript k denotes the shift n& — > + 1. 
However, up to the author's knowledge, no such discretization was known 
till now. The discretization proposed by Schief and Bobenko [5, 6, 7] is 
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very natural from the geometric point of view, but it is given by a three- 
component system rather than one scalar equation: 



AB 

hhi 2 (hih 2 - hi — h 2 ) + h\ 2 + h — 1 j—hih 2 hi 2 = 0, 

A 2 _h\ Bi _ h 2 , , 

A ~ h' B ~ h' [ ' 

The discretization through the permutability properties of the Backlund 
transformations [8, 9, 10] also leads to equations with several components, 
because the Tzitzeica equation does not admit Backlund transformations of 
first order (in contrast to the sine-Gordon equation). 

On the other hand, the discrete theory is usually in perfect parallel with 
the continuous one and therefore one should expect the existence of some 
quad-equation analogous to (1). The aim of this paper is to study the 
following candidate for this role: 

hhi 2 (c~ l h\h 2 — hi — h 2 ) + hi 2 + h — c = (3) 

or, in the exponential form, 

eHl +H 2 -C _ e Jfi _ e H 2 _ e C-H-H 12 _ e -H _ e -H 12 _ 

The parameter c 7^ 0, 00 is essential and cannot be scaled, but the values c, 
— c are equivalent under the change of sign h — > —h and c, c -1 are equivalent 
under the reflection h(ni,n 2 ) — > l/h(—rii,n 2 ). The fixed points c = ±1 are 
special: at these points the equation degenerates into 

hh 12 (ht - l)(h 2 -l) = (h- l)(h 12 - 1) (4) 

which is one of the form of the well known discrete Liouville equation. This 

T\T 2 

is a linearizable equation, namely, the substitution h = maps solutions 

TT 12 

of the linear wave equation t\ 2 — n — r 2 + r = into solutions of (4). The 
general solution of the latter is given therefore explicitly by the cross-ratio 

_ (01 -b)(a- b 2 ) 
(a - b)(ai - b 2 ) 

where a = a(n\) and b = b(n 2 ) are arbitrary functions of one discrete 
variable. At c 7^ ±1, no such explicit formula exists, however we will see 
in Section 3 that the above substitution still makes sense and brings the 
equation to the trilinear /bilinear forms. 

The continuous limit is given by the following substitutions: 

ci— >l + ae 6 , h(ni, n 2 ) i-> 1 + (3e 2 h(x, y), x = em, y = en 2 . 
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In the limit e —¥ 0, the terms up to e 5 in equation (3) vanish identically and 
the coefficients at e 6 recover the equation 

(3 2 (hh xy - h x h y ) = 2(3 3 h 3 - 2a 

which is the Tzitzeica equation (1) up to the change h = e H and obvious 
scalings (certainly, if c = 1 then a = and we obtain the Liouville equation, 
as it should be). 

In the next section we discuss the zero curvature representation for equa- 
tion (3). An associated differential-difference flow is derived in Section 4. It 
serves as a Volterra lattice type discretization of the Sawada-Kotera equa- 
tion [11] and its modification given in Section 5 defines a continuous sym- 
metry of the discrete Tzitzeica equation. The role of such symmetries in the 
theory of quad-equations is well known, see e.g. [12]. 



2 Linear problems 

Integrability of equation (3) is based on the zero curvature representation 
which can be conveniently written as a system of second order difference 
equations. 

Statement 1. The discrete Tzitzeica equation (3) is equivalent to the con- 
sistency conditions of the following equations: 

ipll ~ fJ-ipl = -T7Z — ~T C^ 1 ~ /^) + 7 — -(^12 - f^l), 
hi{n — c) h — c 

4>12+1p = Hip! + 1p 2 ), 
^22-^2 = 7-77 %(^2-W0 + T -(^12-/^2) (5) 

h2{h — c) n — c 

where 

/i = c — (c + 1)A, v = c — (c — 1)A _1 
and A is the spectral parameter. 

The proof is obtained by straightforward computation. Clearly, one can 
substitute ^12 from the second equation into the other two and then the 
system can be rewritten in the matrix form 

^1 = L*, ^ 2 = M$> 

where ^ = (?/>, ipi, ip2) T ■ The consistency condition then takes the standard 
form of the discrete zero curvature representation L 2 M = M\L with 3x3 
matrices. This equation is exactly equivalent to (3). 
The spectral problem for the Tzitzeica equation (1) 

I^XX = V X tp X + Xe~ V 1py, lfj X y = e V 1p, tpyy = A ~ 1 € ~ " lj) x + Vy^y (6) 
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can be easily recovered by the continuous limit. In the matrix form we 
obtain ^ x = U^f?, ^ y = where ^ = (ip, ip x ,ip y ) T , and the zero curvature 
representation U y — V x = [V,U] as the consistency condition. 

Recall that equations (6) are geometrically nothing but the Gauss equa- 
tions for the indefinite affine spheres ip : IR 2 — > M 3 (the definite ones corre- 
spond to the elliptic version of the Tzitzeica equation), see e.g. [13]. 

In order to clarify the geometric meaning of system (5), we compare it 
with the Schief-Bobenko discretization of Gauss equations [5, 6, 7]: 

ipu ~ipi = h \ \ (ipi -ip) + 7 J -— ^(^12 - ^l), 
tpi2 + ip = h(ij>i + ip 2 ), 

1p22 ~ 1p2 = ^ h {tp2 + ^3^(^12 ~ -02)- (7) 

In this case the consistency conditions give exactly system (2). Equations 
(7) describe the so-called discrete indefinite affine spheres, a class of discrete 
surfaces ip : Z 2 — > M 3 characterized by two properties: 

1) ip is a discrete asymptotic net, that is, the points ip, ip\, ipj, 1P2, 
1P2 are coplanar for any (711,712) € 1? (here, k denotes the backward shift 
n k ^-n k - 1); 

2) ip is a discrete affine Lorentz harmonic net, that is, all discrete affine 
normals (defined as vectors 1P12 — ipi ~ ^2 + ip attached to the centers of 
plaquettes \(ipi2 + ipi +fp2 + VO) meet in one point, the origin. 

Property 2) is expressed analytically by second equation (7) and we see 
that it is satisfied for system (5) as well. 

Property 1) is expressed by first and last equations (7), with the co- 
efficients specified by use of the consistency condition. We see that this 
property does not hold for the surface ip (711,712) defined by system (5), but 
it holds for the surface obtained by the gauge transformation 

n 2 ) = fj,^ ni u~ n2 ip(n 1 ,n 2 ). 

Certainly, it is possible to rewrite the system in terms of ip, then property 
1) will be satisfied, but property 2) will be distorted. 

Therefore, one can say that system (5) defines a class of discrete surfaces 
which are, up to a simple coordinate-dependent scaling, both discrete asymp- 
totic nets and discrete affine Lorentz harmonic nets, however in contrast to 
the Schief-Bobenko case these properties are not fulfilled simultaneously. 

To finish this section, we notice that elimination of ip2 and ip\2 from two 
first equations (5) brings to the following ordinary difference equation of 
third order with respect to the shift T\\ 

h u (c - ht) 



uipm + ipu = pb(ip\ + uip), u:-- 



huhih — c 
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(certainly, an analogous equation holds for the shift T 2 ). This linear problem 
is considered in more details in Section 4. 



3 Bilinear equations 

Consider the continuous case first. It is well known (see e.g. [6, 14]) that the 
substitution h = —2(logT) xy maps the algebraic form of Tzitzeica equation 

into the trilinear form 

/ Tyy ^" x yy ^"xxyy \ 

(8) 





( T yy 


T X yy 


Txxyy 


4det 1 


T y 


T~xy 


T~xxy 






T x 


T~xx 



The fact which apparently has not been paid attention before is that a couple 
of simpler bilinear equations can be added, 

3\!~xyTxx T x T~xxy) — TyT xxx TT X xxyi 

3( T xy T yy ~ T y T xyy) = T x T yyy ~ TT X yyy (9) 

which are consistent with (8). These follow from the conservation laws 

which admit integration after the substitution: 

= -6(log t) xx + a(z), = -6(log r) yy + b(y). 

Notice that r-function is defined in (8) up to the multiplication by arbitrary 
functions on x and on y, and this freedom can be fixed by setting a = b = 
without loss of generality. Now substituting r again yields relations 

(log = -6(logr) xa; (logT) X j / , (logr^y = -6(logr)j / j / (logr) X j / 

which are (9). 

In the discrete case the picture is very similar. The substitution h = 

TT 12 

brings (3) 

/i/ii 2 (c _1 /ii/i 2 - hi - h 2 ) + h\2 + h - c = 
to the trilinear form 

C _1 r 2 2Ti 2 Tn + TTi 22 Tn 2 + TiT 2 Tn 22 = TnT 2 ri 22 + TlT 22 Tn 2 + CTTi 2 Tu22 
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which can be rewritten as 




(c - C 1 )TTi2rn22. (10) 



The additional bilinear equations 



T11T12 - CT1T112 = CTT\\Yl - TUIT2, 

TV1T22 - CT 2 Tl22 = CTT1222 - T222T1 (11) 



can be extracted from the multiplicative conservation laws 

«2 h h n (c-hi) vi h h 2 2(c — h 2 ) 



(12) 



u hn hnh\h — c v /i22 h22h 2 h 

These admit integration after the substitution: 

u = a(ni) , v = b(n2) , 

T\\T\ T22T2 

and again, since the r-function is defined up to multiplication by arbitrary 
functions on n\ and 712, hence one can chose integration constants a = b = 1 
without loss of generality. Replacing h in the expressions for u, v yields (11). 
In the case of Bobenko-Schief discretization (2) the substitution for h is 
r 2 r 2 

supplemented by A = a — — , B = b — — with arbitrary a(ni), b(n,2), and 

tth rr 2 2 
this leads to the trilinear equation [5, 6, 7] 

T~22 T122 71122' 

det I r 2 t\2 ni2 | = abrf 2 . (13) 

T T\ Tu 

The bilinear equations are not known for this discretization. Clearly, both 
equations (10) and (13) go to (8) under the corresponding continuous limits. 



4 A difference analog of Sawada— Kotera equation 

In the rest of the paper we change the notation: now we will consider only 
one discrete variable n (identified, say, with m) and it is more convenient 
to reserve subscripts for the order of shift along this variable rather than 
to distinguish between one-step shifts along different variables as before. 
For instance, new /14 is the same as old him and u_3 is the same as uyjy. 
The goal of this section is to derive the differential-difference equation 

u t = u 2 (ii 2 «i - U-1U-2) - u(ui - u-i). (14) 



(3 



Its close relation to the discrete Tzitzeica equation will be revealed in the 
next section, however, this lattice certainly deserves study on its own. It can 
be interpreted as a difference analog of the Sawada-Kotera equation [11] 



U T = U S 



+ 5UU XXX + 5U X U XX + 5U 2 U X 



which appears under the following continuous limit at e — > 0: 



( n t) = - H 

v ' ' 3 9 




/ 4 2e J 



) 



x = en. 



Recall that both flows 



Uf = u{u\ — U-\) and u t » = u 2 (u 2 u\ — u-iU- 2 ) 



are very well known integrable equations: the first one is the celebrated 
Volterra lattice and the second one is the modified Narita-Bogoyavlensky 
lattice (see [15, 16, 17] for a detailed theory and the bibliography). In 
both cases, the continuous limit is the Korteweg-de Vries equation Ut = 
U X xx + 6UU X . However, these lattices are not members of one and the same 
hierarchy, that is the flows df and dt" do not commute. Therefore, one 
should not expect that their linear combination (14) is integrable. But it is, 
and we will prove it by constructing the Lax representation. 
The starting point is the difference spectral problem 



This is a special reduction of a general third order problem, in parallel with 
the continuous case of the Sawada-Kotera equation. We write it in the 
operator form 



where T : n — > n + 1 is the shift operator. The isospectral deformations are 
defined by the Lax equations 



DtiQ-^P) = [A, Q^P] ^ P t = BP- PA, Q t = BQ - QA. 



u t = B(T + u)-{T + u)A, B(T 2 — 1) = (T + u)AT - (uT + 1)A 2 . 



We assume that both A and B are difference operators, that is, Laurent 
polynomials in T. The notation like A 2 is used for the operators with shifted 
coefficients, so that T 2 A = A 2 T 2 . The second equation can be solved as 



Ulp 3 + lp 2 = /J,{lpl + UTp). 



(15) 



Pi/) = /jtQip, P = {uT + l)T 2 , Q = T + u 



These equations are equivalent, for the above P, Q, to the system 



follows: 



A = F(T — T ), B = F\T + u(F — F 3 ) — F 2 T 
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and then the substitution into the first equation yields 

ut = TFu + uFT^ 1 - uF 3 T - T^F^u + F 1 - F 2 + u(F - F 3 )u. (16) 

Notice that this equation admits the reduction F* = F, with respect to the 
usual conjugation of difference operators: T* = T -1 , (FG)* = G*F* and 
f* = f for functions /. Let operator F be of the form 

F = f^T k + ■■■ + fWT + / ( fc '°) + f^T- 1 + ■■■ + £fh~ k , k>0. 

Then collecting coefficients at T k+l , . . . ,T in equation (16) brings to the 
recurrent relations 

u k+1 f?' k) - u/f fc) = 0, 
u*/}*'*- 1 ' - uft^ = ft k) - fi k ' k) + uu k (ft k) - /<*■*>), 

»,./f- 7 l: - U ft' j - 1] = i K3) - fi kJ) + UUj (ft j) - / (fcj) ) 

+ u 3 f^ +1) -uf^\ j = l,...,k-l 
and vanishing of the free term yields the desired equation for u 

u th = 2u(f^ - f™) + u 2 (/ (M) - ft° ] ) + /?' 0) - i m - 

We are interested only in a local evolution, that is, the recurrent relations 
must define all /( fej ) as functions of a finite set of it,. It is easy to see that 
the first equation for can be solved if and only if k is odd, namely 

j?(k,k) _ n _ i . . . Uk _ iUk _ 2i k = 2m + 1 > 0. 
In the simplest case k = 1 we find 

/^ 1 )=u_i, = 1 - «_:,«_! 

and the corresponding flow is exactly (14). So, we arrive to the following 
statement. 

Statement 2. Lattice (14) governs the isospectral deformation of the linear 
problem (15) defined by the equation 

ik = A(i>) = (u^T + 1 - u-xU-2 + U- 2 T- X ){T - T~ 1 )(^). (17) 

The higher flows can be computed analogously. At k = 3 one finds 

/( 3 > 3 ) = mu-t, / (3 ' 2) = w + w-u / (3 ' 1} = 1 -U-i(lV + W-l + W- 2 ), 
j(3,0) _ u _ lU _ 2 ( w _|_ _|_ w _ 2 _|_ w _ 3 ) 
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where w = u(l — u±U-i) and 

U t3 = u(wi(w 3 + W 2 + Wl + W) - W-l(w + W-l + U>_ 2 + W- 3 ) 
—u\{wz + W-l) + U-l(wi + W-3)). 

One can check directly that this is a commuting flow for (14) indeed. It 
should be noted that its right hand side is a sum of three homogeneous 
polynomials: 

Utz = p(3) + p(5) + p(7) 5 deg pM = TOj 

where P^ 3 ^ and P^ define respectively the symmetries for the Volterra and 
the Narita-Bogoyavlensky lattices 

u t ' 3 = P^ = u(ui{u2 + ui + u) — u-i(u + u-i + u_ 2 )), 

U t n = P( 7 ^ = U 2 Ui«2'U3U4 + . . . 

5 Modified lattice 

The role of discrete Miura type transformations for lattice (14) is played by 
the substitutions 

M-: u= h ^- kl \ M+: u= ^=^ 
hihih — c h2fiih — c 

(the first one already appeared in the end of Section 2). 

Statement 3. The substitutions map solutions of the lattice 

^ _ hjc - h) / h(c - hi)(c - h_i)(h 2 hi - fc-ifc-g) _ h ^ + h \ ^ 
hihh-i — c \ (h 2 hih — c){hh-ih— 2 — c) / 

into solutions of lattice (14)- 

Proof. One can easily check that if ^ is a particular solution of linear problem 
(15) corresponding to the value of spectral parameter \x = 1/c then the 
quotient h = (f)/4>i is related to the potential u by the substitution M~ . 
Therefore, the time evolution of h can be found by use of equation (17) for 
<f>t which guarantees that h satisfies some modified lattice equation in the 
form of a conservation law 

(log h) t = (T- l)S(hi,h, h-!,h- 2 ) (19) 

where 

S = —j = \{u-iT + 1 - U-1U-2 + ti_ 2 T~ 1 )(0_i - 4>i) 

[l - ^) + U_2{h„ih_2 - 1) + (1 - «_ 1U _ 2 ) - I) . 



U-l 
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Replacing u and some algebra bring to equation (18). 

The second substitution follows either from the involution h — > h , 
c — > c _1 which preserves equation (18) or from the form invariance of linear 
problem (15) with respect to the change n — > —n. □ 

In order to relate lattice (18) with the discrete Tzitzeica equation (3) 
we notice that u/u = h/h2 and compare this with the first of conservation 
laws (12). It is clear that if the discrete variable n in (18) is identified with 
ni in (3) then the map u — >■ u can be identified with the shift T 2 along the 
second discrete variable ri2- This argument is not quite rigorous, since the 
conservation law is not exactly equivalent to equation (3). Nevertheless, the 
following statement holds true. 

Statement 4. The discrete Tzitzeica equation (3) admits the evolution sym- 
metry (18) along any of the coordinate directions n = n\ or n = that is, 
the differentiation Dt in virtue of lattice (18) is consistent with the discrete 
equation: 

A(Q)|q=o = 0. 

The proof of this identity is straightforward, although rather involved: 
Dt(Q) contains variables h in 12 points {n\ + k, 712), {n\ + k,n2 + 1), k = 
—2, . . . , 3, so that 5 copies of the quad-equation T~ 2 (Q) = • • • = T 2 (Q) = 
are used. 

It is worth noticing that an infinite sequence of conservation laws for 
lattice (14) can be extracted from the single conservation law (19) for the 
modified lattice, by means of the classical trick with the inversion of Miura 
transformation as a power series in the spectral parameter [18]. Indeed, the 
equation u = M~(h) can be solved with respect to h as the formal power 
series 

h = —(l + chM+<?h® + ...) 
u 

with the coefficients recursively found as polynomials in Uj\ 
h^ = ui(l — UU2), 

/j( 2 ) = -u 1 n 2 ((l - u 2 u)(l - U3U1) - u 3 u(l - U4U2)), . . . 
Then the coefficients of the expansion 

]ogh = -log(-u) + ch^ + c 2 (h^ - \{h^) 2 ) +... 

provide densities of the conservation laws Dt(p^) = (T — l)a^: 
p(°)=logti, = u(l — UiU_i), ... 
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6 Concluding remarks 



It should be noted that the discrete Tzitzeica equation (3) falls out of the 
classification scheme based on the notion of 3D-consistency. Indeed, this 
notion is essentially equivalent to the existence of a zero curvature rep- 
resentation with 2x2 matrices [19], while the construction in Section 2 
brings to the 3x3 matrices and it is difficult to expect that the size can 
be reduced. Up to the author's knowledge, this is the first example of a 
quad-equation associated with a third-order spectral problem. According to 
[20] , any 'generic' quad-equation can be included into a consistent triple, but 
the situation remains unclear for equations with degenerate biquadratics. In 
this terminology, equation (3) is degenerate, and this example demonstrates 
that this class of equations can be more complicated than it seems at a first 
glance. The alternative approaches to the integrability of quad-equations 
are developed, e.g. in [21, 22, 23] and one may hope that more examples of 
such kind will be discovered by these methods. 

An open problem is to construct the Backlund transformations for equa- 
tion (3). Probably, the commutativity of these transformations can be for- 
mulated as 3D-consistency of equation (3) with some discrete equations of 
second order in shifts rather than quad-equations. 

The exhaustive classification of integrable Volterra-type lattices ut = 
f(ui,u,u-i) was obtained by Yamilov [15, 16]. However, a little is known 
about the higher order lattices, even in the polynomial case, so that the 
study of lattices like (14) is of interest as well. 
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